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We analyze in the framework of pQCD the properties of the single spin asymmetry in heavy flavor 
production by linearly polarized photons. At leading order, the parallel- perpendicular asymmetry 
in azimuthal distributions of both charm and bottom quark is predicted to be about 20% in a wide 
region of initial energy. Using the soft gluon resummation formalism, we have calculated the next- 
to-leading order corrections to the asymmetry to next-to-leading logarithmic accuracy. It is shown 
that radiative corrections practically do not affect the Born predictions for the azimuthal asymmetry 
at energies of the fixed target experiments. Both leading and next-to-leading order predictions for 
the asymmetry are insensitive to within few percent to theoretical uncertainties in the QCD input 
parameters: ttiq, [ir, (if, Aqcd and in the gluon distribution function. Our analysis shows that 
. nonperturbative corrections to a B-meson azimuthal asymmetry due to the gluon transverse motion 

^ ' in the target are negligible. We conclude that measurements of the single spin asymmetry would 

^ . provide a good test of pQCD applicability to heavy flavor production at fixed target energies. 

Presently, the basic spin- averaged characteristics of heavy flavor hadro-, photo- and electroproduction are known 
, exactly up to the next-to-leading order (NLO) |l]-[|. Two main results of the exact pQCD calculations can be 
formulated as follows. First, the NLO corrections are large; they increase the leading order (LO) predictions for both 
charm and bottom production cross sections approximately by a factor of 2. For this reason, one could expect that 
the higher order corrections as well as the nonperturbative contributions can be essential in these processes, especially 
for the c-quark case. Second, the fixed order predictions are very sensitive to standard uncertainties in the input QCD 
parameters. In fact, the total uncertainties associated with the unknown values of the heavy quark mass, mQ, the 
factorization and renormalization scales, /if and (Ir, ^qcd and the parton distribution functions are so large that 
one can only estimate the order of magnitude of the NLO predictions for total cross sections . For this reason, 
it is very difficult to compare directly, without additional assumptions, the fixed order predictions for spin-averaged 
cross sections with experimental data and thereby to test the pQCD applicability to the heavy quark production. 

Since the spin-averaged characteristics of heavy flavor production are not well defined quantitatively in pQCD it is 
of special interest to study those spin-dependent observables which are stable under variations of input parameters of 
the theory [jlQ(] . In this report we analyze the charm and bottom production by linearly polarized photons, namely 
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I. INTRODUCTION 



X 

the reactions 



1 1 +N ^ Q(Q) + X. (1.1) 



We consider the single spin asymmetry parameter, A(s), which measures the parallel-perpendicular asymmetry in the 
quark azimuthal distribution: 

A( s = 1 da(s,ip = 0) ~ dcr(s,ip = tt/2) 
W V^da(s,(f = 0)+da(s,(p = n/2)' [ ' 1 

Here da(s,(p) = ^-(s,<p) , P 7 is the degree of linear polarization of the incident photon beam, y/s is the centre of 
mass energy of the process ([□]) and ip is the angle between the beam polarization direction and the observed quark 
transverse momentum. 



'Contribution to XV International Seminar on High Energy Physics Problems, Dubna, Russia, Sept. 25-29, 2000 
^E-mail: nikiv@uniphi.yerphi.am 



1 



The properties of the single spin asymmetry at Born level as well as the contributions of nonperturbative effects 
(such as the gluon transverse motion in the target and the heavy quark fragmentation) have been considered in |l0| . 
Using the methods of Refs. |ll-jl4| for the threshold resummation of soft gluons, we have also calculated the next- 



to-leading order corrections to A(s) at next-to-leading logarithmic (NLL) level |lq |. The main results of our analysis 
can be formulated as follows: 



At fixed target energies, the LO predictions for azimuthal asymmetry ( |l.2| ) are not small and can be tested 
experimentally. For instance, 

A(s = 400GeV 2 ) |£° arm « A(s = 400GeV 2 ) |*£ ttom « 0.18. (1.3) 



• Radiative corrections practically do not affect the Born predictions for A(s) at fixed target energies. 

• At energies sufficiently above the production threshold, both leading and next-to-leading order predictions for 
A(s) are insensitive (to within few percent) to uncertainties in the QCD parameters: ttiq, [ir, fip, Aqcd and 
in the gluon distribution function. This implies that theoretical uncertainties in the spin-dependent and spin- 



averaged cross sections (the numerator and denominator of the fraction (1.2), respectively) cancel each other 
with a good accuracy. 

• Nonperturbative corrections to the 6-quark azimuthal asymmetry A(s) due to the gluon transverse motion in 
the target are negligible. Because of the smallness of the c-quark mass, the fcr-kick corrections to A(s) in the 
charm case are larger; they are of order of 20%. 

We conclude that the single spin asymmetry is an observable quantitatively well defined, rapidly convergent in 
pQCD and insensitive to nonperturbative contributions. Measurements of asymmetry parameters would provide a 
good test of the fixed order QCD applicability to heavy flavor production. 

II. SINGLE SPIN ASYMMETRY AT LEADING ORDER 

At leading order, 0(a em as), the only partonic subprocess which is responsible for heavy quark photoproduction is 
the two-body photon-gluon fusion: 

1 (k J )+g(k g )^Q( PQ ) + Q(p-). (2.1) 
The cross section corresponding to the Born diagrams is Jl0| : 

}2~Born P * n ^ n B l\ X X + £ 2 2(1 - /3 2 ) (/3 2 - X 2 ) 



d 2 <7 Born . ~ 2 e 2 Q a em a s (^ B ) 

—— (S, t, if, H R ) = C ; 



l-X 2 {l-X 2 ) 



(1 +7> 7 cos2^) 



(2.2) 



where P 7 is the degree of the photon beam polarization; tp is the angle between the observed quark transverse 
momentum, pq,t, and the beam polarization direction. In (2.2) C is the color factor, C = Tp =Tr(T a T a ) / \N 2 — 1) = 



1/2, and eQ is the quark charge in units of electromagnetic coupling constant. We use the following definition of 
partonic kinematical variables: 

s = (k 1 + k g ) 2 ; i=(k n -p Q ) 2 ; 



u={k g -p Q ) 2 - x = l + 2—^—- (2.3) 



t — m 



where m is the heavy quark mass. 

Unless otherwise stated, the CTEQ5 jl6| parametrization of the gluon distribution function is used. The default 
values of the charm and bottom mass are m c = 1.5 GeV and mi, = 4.75 GcV; A4 = 300 MeV and A5 = 200 MeV. 
The default values of the factorization scale [if chosen for the A(s) asymmetry calculation are \if |charm= 2m c for 
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the case of charm production and /if | Bottom = m b for the bottom case [|9[|i7|]. For the renormalization scale, /ir, we 
use ^ R = n F . 



Let us discuss the pQCD predictions for the asymmetry parameter defined by (1.2). Our calculations of A(s) at 
LO for the c- and &-quark are presented by solid lines in Fig.l. One can see that at energies sufficiently above the 
production threshold the single spin asymmetry A(E 7 ) depends weekly on E~, _E 7 = (s — m^)/2mjv. 
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FIG. 1. Single spin asymmetry, A{E- 1 ), in c- and b- quark production as a function of beam energy Ey = (s — m%) /2mjv; 
the QCD LO predictions with and without the inclusion of fcr-kick effect. 

The most interesting feature of LO predictions for A(E~) is that they are practically insensitive to uncertainties 
in QCD parameters. In particular, changes of the charm quark mass in the interval 1.2 < m c < 1.8 GeV affect the 
quantity A(E^) by less than 6% at energies 40 < E 1 < 1000 GeV. Remember that analogous changes of m c lead to 



variations of total cross sections from a factor of 10 at E~ 



40 GeV to a factor of 3 at £L = 1 TeV. The extreme 



choices mj, = 4.5 and m& = 5 GeV lead to 3% variations of the parameter A{E~) in the case of bottom production 
at energies 250 < E 7 < 1000 GeV. The total cross sections in this case vary from a factor of 3 at E 7 — 250 GeV to a 
factor of 1.5 at E^ = 1 TeV. The changes of A(E 1 ) are less than 3% for choices of /ip in the range \m\, < /j,p < 2mi,. 
For the total cross sections, such changes of hf lead to a factor of 2.7 at E 7 — 250 GeV and of 1.7 at E^ = 1 TeV. We 
have verified also that all the CTEQ3-CTEQ5 versions of the gluon distribution function jl6| as well as the CMKT 
parametrization fisf| lead to asymmetry predictions which coincide with each other with accuracy better than 1.5% 

III. NONPERTURBATIVE CONTRIBUTIONS 

Let us discuss how the pQCD predictions for single spin asymmetry are affected by nonperturbative contributions 
due to the intrinsic transverse motion of the gluon in the target. In our analysis, we use the MNR model (Tt]] 
parametrization of the gluon transverse momentum distribution, 



k g = zk N + k g>T , 

where kjy is the target momentum in the 7./V centre of mass system, 
momentum, k g ^T , has a random Gaussian distribution: 



1 d'N 1 

^Vd^"^f) exp 



(4) 



According to |lj 



(3.1) 

the primordial transverse 
(3.2) 



where k\ 



k 2 



It is evident that the inclusion of this effect results in a dilution of azimuthal asymmetry. In 
p~7| , the parametrization (|]^) (so-called fc^-kick) have been used to describe the single inclusive spectra and the 
correlations. It was found that in charm photoproduction 0.5 < (fcf.) < 2 GeV 2 . 
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Our calculations of the parameter A(s) at LO with the fc^-kick contributions are presented in Fig.l by dashed 
{{14) = 0.5 GeV 2 ) and dotted ((fc 2 ) = 1 GeV 2 ) curves. 

So, we can conclude that nonperturbative corrections to the &-quark asymmetry parameter ( |l.2| ) due to the fey-kick 
effect practically do not affect predictions of the underlying perturbative mechanism: photon-gluon fusion. 

Our calculations of the pr- and x p -distributions of the azimuthal asymmetry in heavy quark photoproduction are 
given in flOI. 



IV. NEXT-TO-LEADING ORDER CORRECTIONS 

The perturbative expansion for the 73 cross section, a{p, /i 2 ), is usually written in terms of scaling functions: 

( 2\ 2 00 k 2 

W) = aemas{ ? )e ^ (4 TOs ( M 2 )) fc ^c^)( P )W^, (4.1) 



k=0 



1=0 



where p = 4m 2 /s, s = (fc 7 + k g ) ; c^ 0,QS> {p) corresponds for the Born contribution (2.2), ^^{p) and c^ 1,0 \p) describe 



the order-ag corrections, factorization scale {p,p = pn. = p) dependent and independent, respectively. 

To calculate the NLO corrections to the single spin asymmetry A{s), we need to take into account the virtual 



0{a em a 2 s ) corrections to the Born process (2.1) and the real gluon emission in the photon-gluon fusion: 



7(M + 9{k g ) Q(PQ) + QiPn) + 9{Pg)- 



(4.2) 



The scale dependent term c^ l,1 \p) which is the coefficient of In ^ can be expressed explicitly in terms of the 

Born cross section using the renormalization group arguments || . Contribution of the scale independent cross secton, 
r (i,o) 



(p), near the threshold can be obtained with help of the soft gluon resummation method J12|— pL4|| . To the next- to 
leading logarithmic accuracy, the soft gluon contribution to the photon-gluon fusion can be written in a factorized 
form as: 



" 2 dSk * BBmn «i) +E {^r) n R{n) Mi, fix)} , 

where 1 1 = t — m 2 , ui = u — m 2 and B Bmn (s, i\, u\) describes the Born level 73 cross section: 



(4.3) 



D Born t - 7 ~ \ 2 
B {S,ti,Ui) = eQOi em OtS 



_L _)_ V± 
Ui ii 



tiUi \ t\U\ 

At NLO, 0{a em ag), the soft gluon corrections to NLL accuracy in the MS factorization scheme are (cf. Ref. fl4|| ) 



(4.4) 



(s,^,^) =2C A 



In (s4/m 2 ) 



.s 4 



C A In 



ReL,g — In 



- 2C F (ReL/3 + r 



5 (s 4 ) C A In [ hn 1 M " 



(4.5) 



where Ca = N c = 3, Cf = 



2N C 



in (4^) are defined as: 



|, §4 = s + ti + u\ and (3 — y/1 — p. The function Lp and the plus-distribution 
1 - 2m 2 /s I", (\ - s 



Lf3 = 







In 



1 + 



(4.6) 



V {s A /m 2 )~ 


= lim < 


h 





In' {si/i 



Si 



' (s 4 - A) 



l + l 



111 



1+1 



A 



S{s±) 



(4.7) 
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Eq.(4.5) describes very well the threshold behavior of the exact partonic cross sections c^'°\p) and c^ 1,:L ^(p) and 
gives, at NLO, practically whole contribution to the unpolarized bottom production in photon-hadron reactions at 
fixed target energies, E 1 < 1 TeV. 

All terms in the r.h.s. of Eq.( |4.5| ) originate from the so-called collinear p g x — ► 0, and soft, p g — ► , components of 
cross section. Since the azimuthal angle ip is the same for both 75 and QQ centre of mass systems in the collinear and 



soft limits, Eq.(4.3) can be generalized to the spin-dependent case substituting the spin-averaged Born cross section 
by the (^-dependent one: B Boin (s,ti,ui) — > B BoIn (s,ii,Ui,(p). 

The results of our calculations of the single spin asymmetry A(s) at NLO to NLL accuracy in c- and 6-quark 
production are presented by dashed line in Fig. 2. The details of calculations as well as the higher order predictions 
for A(s) will be given in |i~5f|. 
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FIG. 2. Single spin asymmetry, A(E y ), in c- and &-quark production at LO (solid curve) and at NLO to NLL accuracy 
(dashed curve). 



One can see from Fig. 2 that the NLO NLL and Born predictions for A(s) coincide with each other with accuracy 
better than 2%. We have verified that the azimuthal asymmetry is independent (to within few percent) of theoretical 
uncertainties in the QCD input parameters (ttiq, p,R, p,p and Aqcd) at NLO too. 



V. CONCLUSION 



Our analysis shows that the NLO corrections practically do not affect the Born predictions for the single spin asym- 
metry in heavy flavor production by linearly polarized photons at fixed target energies. So, the quantity A(s) is an 
observable quantitatively well defined, rapidly convergent and insensitive to nonperturbative contributions. Measure- 
ments of the azimuthal asymmetry in bottom production would be a good test of the conventional parton model 
based on pQCD. 
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